A warning on the determination of the halo mass 



Darfo Nunez, Alma X. Gonzalez-Morales*, Jorge L. Cervantes-Cota^** 

and Tonatiuh Matos*'® 

* Institute de Ciencias Nucleares, Universidad National Autonoma de Mexico, A.P. 70-543, 04510 

Mexico D.F., Mexico 
^Depto. de Fi'sica, Institute National de Investigaciones Nucleares, Mexico, 
** Berkeley Center for Cosmological Physics, University of California, Berkeley, CA, US 
^Departamento de Fi'sica, Centro de Investigation y de E studio s Avanzados del IPN, A.P. 14-740, 

07000 Mexico D. F, Mexico 
^Institute Avanzado de Cosmologi'a, IAC 

Abstract. We summarize our studies on the determination of the mass of the dark matter halo, 
based on observations of rotation curves of test particles or of the gravitational lensing. As we show, 
it is not uncommon that some studies on the nature of dark matter include extra assumptions, some 
even on the very nature of the dark matter, what we want to determine!, and that bias the studies 
and the results obtained from the observation and, in some cases, imply an inconsistent system 
altogether. 
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INTRODUCTION 

Einstein's equations (including the cosmological constant) have proven to accurately 
describe the Universe at the Solar System scale (Precession, GPS, and gravitational 
lensing, the general relativistic experiment per excellence), as well as at large scales, 
100 Mpc, where the Friedmann-Robertson-Walker spacetime determines the standard 
cosmological model, which is observationally corrobotated by the Cosmic Microwave 
Background Radiation; by the redshits in galaxies, and by the relative abundance of 
Hidrogen to Helio, which is 1 : 4. 

Based on these facts, our position regarding the new set of observations, is to consider 
the Einstein's equations to be valid at all the intermediate scales, without any modifica- 
tion, and try to understand its consequences. 

These new set of cosmological observations are the following four main groups: 
the acoustic peaks in the cosmic microwave background radiation, Supernovae type la 
data, rotation curves of spirals and dynamics of galactic clusters, and their gravitational 
lensing [1, 2, 3]. The last two of these observations are related most significantly with the 
presence of dark matter (DM), where a consistent model considers a Dark Matter halo 
surrounding the galaxies and galactic clusters. All of these cosmological observations 
are consistently described by the ACDM model. We mean by this that, as we consider 
that our hypothesis is that the Einstein's equations do describe the dynamics of the bodies 
moving on a curved background, which in turn is curved due to the presence of matter, 
then, when a given observation, say the rotational curve profile, is not explained by the 



amount and disctribution of matter that we see, there must be extra matter that we are 
not seeing, but that affects the motion as the Einstein's equation dictates. 

There has been many attempts to model DM halos, some of them have shown that 
general relativistic effects can be important [4, 5, 6, 7, 8, 9, 10, 11], but in general, it is 
fair to say that the make use of approximations or limits that in some cases represents 
assumptions about DM properties. Even lensing images or distortions of background 
galaxies due to the space-time curvature [12, 13] are usually described within the 
Newtonian framework, which is remarkable as long as lensing is a purely relativistic 
effect. The point is that usually assumptions are made precisely on the nature of dark 
matter, which is what one is trying to determine, and those extra assumptions can bais 
the conclusions of produce an altogther inconsisten description. 

The nature of the Dark matter is unknown. The scalar field DM model is an alternative 
proposed in the past [14? ] to fit the observed amount of substructure [15], the critical 
mass of galaxies [16], the rotation curves of galaxies [17], the central density profile of 
LSB galaxies [18], the evolution of the cosmological densities [19], among other topics. 

In the present work we focus on a given simple model for the DM halo and explicitly 
show what exactly is determined by the observations and what comes as extra assump- 
tions. This is particularly important as the unknown nature of the DM is one of the 
most relevant questions that one would like to solve. It is clear that to make assumptions 
specifically on the nature of DM, in order to obtain information on the nature of the DM, 
is skewing the problem. 



THE MODEL 

In our study we consider a static and spherically symmetric space-time in General 
Relativity, described by the line element: 

ds 2 = _ e 2*/c* c 2 df 2 Ijtf (1) 

c 2 r 

where dQ. = dO 2 + sin 2 dd(p 2 . The gravitational potential 4>(r) and the mass function 
m(r) are functions of the radial coordinate only. In fact, due to the symmetries of this 
space-time, all physical quantities depend only on r. The Einstein's equations are the 
known set [20]: 

m' 
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where prime ' = d/dr. The potentials of such space-time are determined by the DM 
halo. In order to see the difference between two hypothesis on the nature of DM, we will 
consider two types of composition for halos: a perfect fluid and a scalar field. 

The above equations are complemented by the conservation equation of the matter- 
energy generating the curvature of the space-time. But given the different nature of the 
fluids considered here, this equation is treated separately for each fluid. 
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Perfect fluid 



In the case of the perfect fluid, the stress-energy tensor is given by v = (pc 2 + 
p)u^u v + pg^, v , where the density is p = (1 +e)p , where Po is the rest mass energy 
density and e the internal energy per unit mass, is the co-movil four velocity, 
normalized as = — c 2 , and p is the pressure. The conservation equation, T^ v;!1 = 0, 
implies the field equation: 

(pc 2 + p)* +p' = 0, (4) 

which can be rewritten as 

T r r'+(1?-If)*=0. (5) 



Scalar field 

Now, for the scalar field the stress energy tensor is given by 

r MV = Mv-\§nv (g af *M,t+iv{*)) , (6) 

where ;Ce = d§ /dx a ; and V (0) is the scalar potential. The components of the stress- 
energy tensor are 

From the conservation equation for the scalar field, r^ v;jU = 0, one obtains a field 
equation, the Klein-Gordon equation, 



(8) 



that can be written as 
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(9) 




which is remarkable similar to the field equation for the perfect fluid, Eq.(5). Given this 
similarity, it is convenient for our mathematical description to consider the single field 
equation for both types of matter 
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in which a = for the perfect fluid, and a = 1 for the scalar field. Notice that if one 
considers a sort of perfect fluid given by T^ v = diag(—p,p,pi,pi),pi (some times called 
"tangential" pressure) denotes a term representing the ignorance we have on the features 
of the fluid. This presure pi is related to the other fluid variables as pi = (1 — a) p — ap, 
(see Eq. (10)), where "a" takes, in principle, any value. There are works which have 
discussed this field equation considering a as a free parameter [11, 9, 10]. For the 
purpose of the present work we will consider only the two extremal cases, a = and 
a= 1, but the discussion can be directly applied for these cases as well. 

In this way, the system of equations which must be solved are, the Einstein's equa- 
tions, Eqs. (2,3), and the field equation, Eq. (10). In either case, there are four unknown 
functions, m, 4>, p and p, for the case where the curvature of the space-time is due to the 
perfect fluid, or m, 4>, (j) and V (0) when the curvature is caused by the scalar field. Thus, 
we have three equations for four unknown functions. In either case, we need only one 
extra data. It is important to underline this fact. Once the extra data is given, there is no 
more room left for any other assumption, the rest of the functions are determined by the 
system of equations. If, for instance, we give an equation of state for the perfect fluid, 
p = p(p) or, in the case of the scalar field, an explicit form for the potential, V (0), there 
is no freedom left to choose the form of the rest of the functions, they will be determined 
by the system of equations. 

Following the line of work presented in [6], we use observational results to close 
the system of equations. In the case of galactic halos, two main observations can serve 
to obtain the desired information: measurements of rotation curves in spirals and light 
deflection by lensing. In this work we choose the former to complement the above field 
equations and use the latter to discriminate between different halo type models. 



ROTATION CURVES 

The motion of test particles in such spacetime is determined by the geodesic equations 
and, for test particles in circular motion, there is a relationship between the gravitational 
potential, 3>, and the tangential velocity of those particles, v c : 

4>' B 2 
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where we have defined /3 2 = This tangential velocity is the one measured by 
observations of rotation curves in galaxies. Thus, v c is an observable function, and 
by means of Eq. (11), the gravitational potential can be determined. Thus, given this 
observable, there is no room left for an equation of state for the perfect fluid or for a 
given scalar field potential. 

Moreover, as long as the magnitude of the observed velocities are small with respect 
to the speed of light, this justifies the validity of one of the weak field approximations 
4>/c 2 << 1 that one usually assumes by taking the weak field limit. Here we want to 
emphasize that the approximations 2Gm/c 2 r « 1 and especially p « p are, in 
general, extra hypothesis which strongly depend upon the nature of the DM type. It 
is clear that if all these conditions are satisfied, then the above system of equations, Eqs. 



(2, 3, 4), reduces to the hydrodynamic set of equations for the case of the perfect fluid 
model. But, for example in the case of the scalar field, there is no Newtonian limit, and 
one has to be careful with these approximations. 

After substituting Eq. (11) into the gravity equations, Eqs. (2, 3), we obtain an 
equation (with no approximations) for the mass function as the only free function 
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m' + P(r)m = Q(r) (12) 
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The functions P(r) and Q(r) depend on the type of fluid we are dealing with (a) and on 
the rotation curves profile. 

The mass function can be expressed in terms of the gravitational potential, <£, through 
the integral 

jernr'w Q{r)dr+c 

e f r P(r')dr< • KL * } 

where the value of the integration constant, C, is set by the appropriate boundary 
conditions. 

For the case of the perfect fluid, the density and pressure are directly computed from 
Eqs. (2) and (3), respectively. 

For the scalar field, using the expressions Eqs. (7), we obtain that 

4 / Gm' Gm \ 

"(♦w) = 4?^ m{l+2 ^ +m '^-ik{^)' <16) 

Once the function /3 (r) is given, the scalar field and its potential are straightforwardly 
determined in terms of the radial coordinate. In order to obtain the form V (0), one needs 
to invert the solution for the scalar field (r = r(0)), and to substitute it into Eq. (16). As 
shown below, this procedure works at least for simple /3 functions. 

In this way, we have shown that the mass function m, associated to a galactic halo 
by means of the rotation velocity strongly depends on the DM model that is being 
considered. The single observation of the rotation curve is not sufficient to determine the 
nature of DM and hence the mass associated with the halo. Moreover, we have shown 
that the relationship between the pressure and the density, or between the scalar field 
and the scalar potential is fixed, up to integration constants, once the rotation velocity is 
employed. 



LENSING 



The other observation concerns the gravitational lensing, that for the line element given 
by Eq. (1), the deflection of the light ray, A(p, at the radius of maximal approach, r m , is 
given by [21], 



Since the gravitational potentials and the fluid variables are already determined by 
the rotation curves of spirals, deflection angle measurements can serve to discriminate 
between models. Here we deal with two examples, perfect fluid and scalar field DM 
models. Yet, observations of spirals that lenses light are not very common, however the 
first examples of them has recently appeared [22]. We remark that from the expression 
of the deflection angle, it is a large step to infer the mass of the DM halo based solely 
on the observation of the deflection angle. A supposition has to be made on the relation 
between the gravitational potential, 4>, and the mass function, m [21]. Such supposition, 
as we have shown, not only strongly depends on the type of matter, but also on the 
specific characteristics of the type of matter considered. 

In the next section we present some examples of known rotation curve profiles to give 
a quantitative description to these conclusions. 



The idea in this section is to stress the conclusions that we are presenting by means 
of considering a typical observation of rotation velocities in spirals and to directly 
determine the gravitational mass in each case, when the DM is a perfect fluid (dust) 
and when it is a scalar field. 

In practice we can consider a velocity distribution, as a phenomenological model, for 
instance the velocity profile coming from N-body simulations given by NFW [23, 24] 
or a Burkert profile [25] given by the phenomenological of rotation curves [26], to 
determine the mass of each type of fluid. We will show that the gravitational mass 
inferred by the same velocity profile is strikingly different for the perfect fluid and scalar 
field cases. 



We will consider the simplest case of constant rotation curves as our first example. 
Although there are some examples of galaxies that present a constant velocity profile, for 
a few disk length scales [27], this is not a typical behavior, being our own Galaxy a good 
counter example [28] and, in fact, there is an important rotation curve phenomenology 
described by the Universal Rotation Curve [26, 29]. However, the constant velocity 
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EXAMPLES 



Constant velocity profile 



profile offers us the mathematical simplicity to obtain straightforward analytical results 
and to show the main point of our work. 

For the gravitational potential, from Eq. (11), when the velocity function is a constant, 
A), we get 

/ r \ ft 2 

<D = c 2 lnf — J . (18) 
The mass function can be analytically obtained for any value of the parameter a as: 

c 2 j /3() 2 (2(1 -fl) -ft) 2 ) (l+2^)(3-2a-^y 

m Po = - I — ^ —> r ±cr ^«+A> 2 I , (19) 
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where C is the integration constant of Eq. (14). For the case of DM described by a perfect 
fluid, a = 0, we fix this constant to zero in order to avoid changes in the signature of 
the line element, Eq. (1). Thus, the mass function, and the corresponding pressure and 
density in the case of the perfect fluid are given by 
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= const. (23) 



We can see in the limit of very small velocities, /3o << 1, we recover the Newtonian 
limit, and pressure is negligible with respect to the density as we mentioned above, 
however, it is not zero and, actually, we obtain a barotropic equation of state, p = wop. 

On the other hand, considering the DM halo due to a scalar field, the mass function is 
obtained from Eq. (19), with a = 1. In this case, the mass function has a very peculiar 
behavior. The first term is small, proportional to A) 4 , but negative. The second term, 

proportional to the constant C, goes as r~( 1+2 ^° ), thus, by choosing a positive value 
for the constant C, one can have a positive mass function for a large region, but this 
function will present a divergence at the origin. Of course this result was expected, 
as the space-time metric is static. In order to avoid this problem, we had to take non- 
static space-times, like the oscillatons [16], but this is beyond the scope of this work. It 
can be shown, however, that this divergence is covered by an apparent horizon. Some 
features of this case of scalar field with a non zero constant C in the mass function, have 
been discussed in [9]. For the purpose of this work, we only notice that the geometric 



functions and those of the scalar field, have a non intuitive behavior, but are consistent 
with the rotational curve. Explicitly, for the case of C = 0, the mass function is 



m?/ = "^T^ r ' (24) 

and, with the geometric functions determined, the scalar field and scalar potential are 
completely fixed, given by: 



C 4 2/3n 2 =c9 

^(0) = e T2 V^^. (27) 



^ G (l-^0 2 )ro 2 



where the expression for the scalar field, 0(r), was inverted to obtain r(0), and then 
express the scalar potential in terms of 0, as explained previously. The "effective mass" 
of the scalar field, m e ff ~ , 2 2 , depends inversely on the characteristic distance of 

the halo. This distance is of the order of kilo-parsecs, and /3o ~ 10~ 3 , thus it will turn into 
a typical mass for scalar field in a galaxy, which corresponds a very light boson mass 
~ 10~ 23 eV/c 2 . This result is in agreement with the one obtained in previous works, see 
for example [15]. 

Going back to our previous discussion, notice how remarkably different are the mass 
expressions derived from each type of fluid, being both consistent with the observed 
rotation velocities. This is the simplest case in which we can show how the single 
observation of the rotation velocities in halos determines the features of the perfect fluid 
model or the scalar field. 

Although the mass associated to the scalar field results negative and this can be taken 
as a no-go result for static scalar field halos [30], rotation curves of spirals are not 
exactly flat (see discussion in [26, 29]) and, in addition, we have to be cautious with 
the supposition of a static metric which is very restrictive for the scalar field . Thus, the 
above-result should not be taken as definitive, at most, it should be taken as a remark 
that a static DM halo is not well described by a static scalar field. A negative mass, 
or positive gravitational potential, is known since long time ago [31] from the fits to 
rotation curves using modifications of newtonian gravity in which a scalar field induces 
a Yukawa-type force. At the end, demanding a constant velocity profile all the way in 
the radial direction implies an effective repulsive force to be acting on test particles in 
the galaxy. 

In any case, it emphasizes our point in showing how strongly depends the determina- 
tion of the mass of the DM halo on the type of matter considered to describe it. 



The deflection angle for the case of constant rotation velocity, considering the perfect 
fluid and the scalar field with the constant C = 0, implies the following expressions: 
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where we have defined x = ^-f. Since in any case the deflection angle is a constant, i.e. 
it does not depend of the maximal approach radio r m , it can be evaluated for a typical 
value of the velocity. For comparison we take the value /3o = 1 /1200, that corresponds 
to a velocity of v c = 250 km/s. Evaluating the deflection angle, we get 



for both cases the deflection angle is given in arc second units. We can see that there is a 
difference of almost half arc second between them, and the simultaneous observation of 
the rotation velocity and the deflection of light produced by the galactic halo, can teach 
us about the true nature of the DM. 

We notice that the deflection angle for the scalar field with a non-zero value of the 
constant C in the mass function takes very large values, a fact which certainly allows us 
to discard this option as a model for the DM halo, independently of any interpretations 
of the mass function. 

Now we study an example that is less striking though. 



Independently from its origin, the NFW profile [23, 24] is considered per se as a 
viable fitting model to describe galactic kinematics. This profile has been subject to 
geometrical studies elsewhere [32, 33, 32]. In this example, we assume this profile 
as a valid phenomenological galactic profile for the galactic data. We obtain the usual 
expression for the mass derived within this description, and compare it with the same 
form of the rotation velocity, but considering that it is due to a DM halo composed of a 
scalar field. 

The velocity profile in the NFW model [23, 24] is given by 



cc p f = 0.899547, 
a sf = 0.449546, 



(31) 
(32) 



Navarro-Frenk- White (NFW) velocity profile 




(33) 



where Oq = AkGPqtq 2 is a characteristic velocity of stars in the halo, given in terms of 
a characteristic density, and r$ is a scale radius. Given this velocity profile, we have to 



solve Eq. (14) with a = for the perfect fluid and with a = 1 for the scalar field. In neither 
case there is an analytical solution, thus we have integrated the equations numerically. 
We do not want to treat here specific galaxies but to emphasize the differences between 
the galaxy models. Therefore, we set Co and tq to some typical values. In our plots 
we assume geometric units (G = c = 1), and therefore the characteristic velocity takes 
values, < Co < 1, and the mass is less than the unity. For definiteness, we assume 
Co = 0.001 and tq = 1. In figure (1) we plot both halo masses (perfect fluid and scalar 
field). Disregarding the behavior near the origin, as long as we are considering the 
outside region, as mentioned above, we see that the mass associated to the halo in each 
case are different. We now consider lensing. By integrating Eq. (1 1) for the given rotation 
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FIGURE 1. Comparison of the masses using the rotation velocity profile from the NFW model with a 
perfect fluid (upper curve) and scalar field (lower curve). 

velocity, Eq. (33), we obtain the following expression for the gravitational potential: 

lnfl + f) 

<J> = -c r ^ ^. (34) 

r 

We substitute this expression, together with the corresponding numerical solution one 
for the mass in each case, in the equation for the deflection angle, Eq. (17), and perform 
the integration varying the value of the radius of maximal approach, r m . The results are 
plotted in figure 2. As we see, the observation of the deflection angle can determine 
which type of matter is actually composing the DM halo. 
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FIGURE 2. Deflection angle generated by the gravitational lensing of a NFW rotation profile with a 
perfect fluid (upper curve) and scalar field (lower curve). 

In this example, the mass associated to the scalar field model is essentially positive 
and a well-behaved function, that is, it does not follow the no-go result mentioned in 





the previous example. Though there is a small region where the mass becomes negative. 
This is due to the fact that we are demanding the velocity profile to grow in that region. 
In a real setting however the stellar disc adds to the velocity profile, thus we expect that 
its contribution avoids negative mass regions for the scalar field. 

With these examples it is clear how two different types of matter (perfect fluid and 
scalar field) can be consistent with the observation of rotation curves of DM halos, 
though they lead to different conclusions to the mass function inferred. The deflection 
of light can then be used to discriminate between the two models. Even though the mass 
function for some model has not an intuitively expected behavior, it is necessary to use 
the observation in order to discard the model, being aware of the assumptions made 
during the derivation of such conclusions. 
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